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SCHOOL  OF  CIVIL  AND  ENVIRONMENTAL  ENGINEERING 
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June  30,  1997  • 

Capt.  Michael  Chipley 
Air  Force  Office  Of  Scientific  Research 
Dir^torate  of  Aerospace  &  Materials  Science 
1 1 0  Duncan  Avenue,  Suite  B 1 1 5 
BoUing  AFB,  DC  20332-8080 


Dear  Mike: 

Enclosed  ple^e  find  the  final  report  for  our  AFOSR-sponsored  research  project  entitled 
“Vapor  Transport  in  Unsaturated  Soil”  (project  no;  F49620-94-1-0157).  As  we  discussed  in  a 
recent  exchange  of  e-mail  messages,  several  aspects  of  the  work  on  this  project  are  ongoing  under 
the  auspices  of  our  augmentation  grant  “Augmentation  for  Basic  Research  on  Vapor  Transport  in 
Soil”  (project  no:  F49620-94-1-0316).  As  a  result,  we  have  not  completed  the  full  analysis  of  our 
data  and  its  implications.  The  aspects  of  this  research  that  have  been  wrapped  up  are  fully 
described  in  the  Cornell  University  M.S.  thesis  of  Mr.  Sean  Lorden,  and  representative  portions  of 
that  material  have  been  extracted  for  this  report  A  publication  based  on  Sean’s  research  is 
currently  in  preparation.  Three  additional  graduate  students  are  currendy  involved  in  the  research. 
Two,  Mr.  Jason  Gilbert  and  Mr.  Brent  Alspach  are  expected  to  complete  their  M.S.  thesis  this 
summer,  and  publications  based  on  their  work  are  anticipated  to  be  forthcoming.  Preliminary 
results  from  their  research  is  summarized  in  this  report.  The  third  student,  Mr.  Kevin  Saunders, 


will  be  continuing  his  research, through  till  the  end  of  the  augmentation  award  (6/30/98). 

By  the  end  of  the  augmentation  award,  the  AFOSR-funded  research  will  have  supported 
the  training  of  four  M.S.-level  environmental  engineers  whose  future  careers  will  contribute  to  the 
solution  of  this  nations  environmental  problems.  In  addition  to  the  M.S.  thesis  of  each  student, 
my  expectation  in  that  much  of  this  work  wUl  be  suitable  for  publication  in  peer  reviewed  journals. 

I  take  this  opportunity  to  thank  you  and  the  Air  Force  Office  of  Scientific  Research  for  your 
support  of  these  investigations;  Basic  research  of  the  type  done  here  serves  as  the  foundation  for 
technology,  such  as  that  used  for  remediation  of  contaminated  soils,  but  also  for  future  applications 
that  we  can  not  yet  anticipate.  The  Air  Force’s  vision  in  supporting  activities  of  this  type  is  to  be 
commended. 

If  any  additional  information  is  needed  regarding  our  research,  or  any  alterations  are  needed 
to  this  report,  please  contact  me  by  phone  (607-255-7571),  fax  (607-255-9004),  or  e-mail 
(LWL3@comell.edu). 

Sincere^^^^ 


Leonard  W.  Lion 
Professor 
enclosure 
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VAPOR  TRANSPORT  IN  UNSATURATED  SOIL 

Abstract 

The  transport  of  volatile  organic  compounds  (VOCs)  is  influenced  by  the  form  in  which  VOCs 
exist  and  the  physical/chemical  properties  of  the  porous  medium.  VOC  transport  mechanisms  in 
the  unsaturated  zone  include  diffusion  and  advection  in  the  gas  phase.  VOC  advection  requires 
mass  transfer  of  sorbed  VOCs  to  the  mobile  gas  phase,  and  given  the  multiple  processes  that 
control  VOC  sorption,  transfer  rates  may  be  best  described  using  probability  distributions.  When 
analyzing  the  retarded  transport  of  VOCs,  it  is  also  important  to  consider  that  isotherms  for  the 
sorption  of  gases  are  typically  non  linear  at  high  relative  pressures,  and  that  VOCs  commonly 
occur  at  waste  sites  as  mixtures  that  can  display  competitive  sorption  effects.  In  order  to  address 
these  issues,  experiments  have.been  performed  to  determine  the  desorption  kinetics,  diffusion,  and 
advection  of  two  common  subsurface  VOC  contaminants:  trichloroethylene  (TCE),  and  chloroform 
(CF)  both  individually  and  as  gas  mixtures. 

Effective  diffusion  coefficients  of  TCE  have  been  determined  in  the  absence  and  presence  of  CF. 
At  high  relative  pressure  the  “effective”  distribution  coefficient  for  these  vapors  decreased  and 
consequently  their  diffusive  transport  was  enhanced.  This  result  is  attributed  to  the  presence  of 
non-linear  sorption  isotherms  for  the  vapors.  Advective  transport  of  TCE  and  CF  also  displayed 
mild  gas  concentration  effects  suggesting  that  the  accurate  prediction  of  vapor  transport  may 
require  consideration  of  reactions  at  the  gas  water  interface  including  vapor  condensation  and 
competition. 

The  breakthrough  and  subsequent  desorptive  release  of  TCE  vapor  from  a  column  of  aquifer  sand 
was  well  predicted  by  a  transport  model  that  employed  a  Gamma  distribution  of  mass  transfer  rate 
constants.  Comparisons  of  experimental  data  and  model  predictions  have  been  successfully  made 
at  two  different  gas  flow  velocities  for  sand  exposed  to  low  (P/Pq  =  10%)  and  high  (P/Pq  =  90%) 
relative  pressures  of  TCE. 
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VAPOR  TRANSPORT  IN  UNSATURATED  SOIL 


Objective; 

The  main  objective  of  this  research  was  to  investigate  the  parameters  that  govern  advective  and 
diffusive  transport  of  volatile  organic  compounds  (VOCs)  in  unsaturated  porous  media,  with 
special  attention  to  vapor  condensation,  interaction  of  vapor  mixtures,  and  mass  transfer  kinetics. 

Background/Introduction: 

As  VOCs  travel  through  a  soil,  they  can  be  bound  by  a  variety  of  mechanisms.  For  example,  it  is 
widely  known  that  nonionic  hydrophobic  compounds  such  as  VOCs  can  partition  from  the 
aqueous  phase  into  soil  organic  matter.  In  the  unsaturated  zone,  water  exists  as  a  film  coating 
mineral  surfaces  as  well  as  bulk  water  condensed  in  cracks,  fissures  and  pores,  and  VOC  uptake 
can  occur  as  gaseous  VOCs  dissolve  into  this  stationary  aqueous  phase.  Additional  mechanisms  of 
VOC  sorption  include  condensation  in  pore  space  not  occupied  by  water  and  adsorption  at  the  gas- 
water  interface.  The  various  VOC  uptake  processes  are  illustrated  in  Figure  1. 

Unlike  partitioning  into  soil  organic  matter  and  dissolution  into  pore  water,  the  adsorption  and 
condensation  reactions  of  VOCs  at  the  gas  water  interface  are  not  necessarily  governed  by  a  linear 
isotherm.  The  adsorption  of  individual  vapors  at  the  gas  water  interface  is  also  expected  to  be 
dependent  on  the  composition  of  the  gas  phase  as  different  VOCs  may  display  differing  affinities 
for  the  surface.  Consequently,  the  magnitude  of  the  sorptive  retardation  of  VOCs  is  anticipated  to 
depend  on  the  concentration  and  composition  of  the  vapor  phase.  Any  variability  in  vapor 
retardation  will,  in  turn,  be  manifested  in  the  transient  advective  and  diffusive  transport  of  VOCs. 


Sorbed  VOC  in  mineral-phase 
cracks,  fissures,  and  pores  ^ 

Dissolved  VOC  in 

sorbent-bound 

moisture  ^  . 

Sorbed  VOC  in 

microporous 

organic  matter 


Condensed 

VOC 


Sorbed  VOC  at 
the  gas-water 
interface 


Figure  1.  Mechanisms  of  VOC  Uptake  on  Unsaturated  Porous  Media. 


Any  mechanistically  based  mathemaUcal  description  of  VOC  transport  in  the  vadose  zone  must 
incorporate  mass  transfer  of  VOCs  from  the  porous  matrix  into  the  mobile  gas  phase.  Each  of  the 
VOC  uptake  mechanisms  indicated  in  Figure  1  will  be  governed  by  different  dynamic  controls. 
The  distribution  of  pore  sizes  and  geometries  in  any  given  porous  medium  will  result  in  a  variety  of 
diffusional  paths  and  diffusive  release  rates  for  VOCs.  Also,  many  different  types  of  organic 
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matter  with  different  VOC  sorption  affinities  and  rates  are  present  in  soils.  Furthermore,  some 
types  of  VOC  uptake  (such  as  condensation  and  adsorption  at  the  gas-water  interface)  are  not  linear 
partitioning  processes,  leading  to  the  expectation  that  there  will  be  nonlinear  isotherms  and 
dynamic  competition  between  VOCs.  As  a  result  of  these  complexities,  an  array  of  VOC  release 
rates  is  expected,  and  the  mechanistic  processes  controlling  sorption/desorption  kinetics  may 
change  over  time.  This  type  of  problem  lends  itself  to  treatment  with  statistical  distributions.  A 
Gamma  density  function  has  been  employed  to  describe  variable  release  rates  of  sorbed 
hydrophobic  organic  contaminants  from  saturated  soil  systems  (Connaughton  et  al,  1994;  Chen, 
1994;  Chen  et  al,  1997;  Culver  et  al,  1997)  and,  given  its  prior  success,  this  model  warrants 
evaluation  for  its  suitability  in  the  modeling  of  VOC  transport. 


Materials/Methods 

A  fme-grained  aquifer  sand  obtained  from  a  quarry  in  Newfield,  NY  was  used  as  the  sorbent  in 
this  investigation.  This  sand  has  previously  been  characterized  with  respect  to  specific  surface  area 
(3.3  m^/g  by  B.E.T.  analysis  with  N2),  organic  carbon  content  (0.91%),  and  particle  size 
distribution  (93.7%  sand,  6.3%  silt)  (Chen  et  al,  1995).  Three  VOCs  were  chosen  for  use: 
trichloroethylene  (TCE),  chloroform  (CF)  and  perchloroethylene  (PCE).  These  VOCs  are 
common  components  of  solvents  and  degreasers,  and  are  suspect  carcinogens.  TCE  is  one  of  the 
most  frequently  detected  subsurface  contaminants  at  hazardous  waste  sites  in  this  country.  Also, 
results  from  preliminary  sorption  studies  conducted  by  Thoms  and  Lion  (1992)  showed 
competitive  sorptive  interactions  between  CF,  TCE  and  PCE.  The  content  of  this  report  is 
principally  confined  to  the  experiments  performed  with  TCE  and  CF. 

In  support  of  the  VOC  transport  experiments  sorption  of  trichloroethylene  (TCE)  was  characterized 
over  the  full  range  of  relative  pressures  onto  the  aquifer  sand  sorbent  at  90%  relative  humidity. 
Non-linear  adsorption  isotherms  were  developed  for  water  and  TCE  using  the  gravimetric 
apparatus  described  by  Ong  and  Lion  (1991).  VOC  sorption  isotherms  were  obtained  at  low 
relative  pressure  under  water-saturated  and  unsaturated  conditions  using  the  method  of  Garbarini 
and  Lion  (1985)  as  modified  by  Peterson  et  al  (1988). 

Diffusion  Experiments 

A  number  of  different  column  configurations  are  commonly  used  in  diffusion  research,  relying  on 
different  modes  of  contaminant  introduction  and  concentration  monitoring  (Batterman  et  al, 
1996;  Mutter  et  al,  1992;  McCarthy  and  Johnson,  1995;  Petersen  et  al,  1994;  Voudrias  and  Li, 
1993).  In  this  research,  a  modified  reservoir-type  system  was  used.  The  diffusion  column  setup 
was  modeled  after  a  system  successfully  implemented  in  previous  work  at  Cornell  (Peterson  et  al, 
1988).  Typically,  in  prior  research  reservoir  setups  have.relied  on  vapor  source  concentrations  that 
correspond  to  the  saturation  vapor  pressure  (Psat)  of  an  individual  compound.  The  apparatus  used 
in  this  study,  however,  allowed  for  variable  VOC  vapor  concentrations  and  for  mixtures  of  VOCs. 
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As  shown  in  Figure  2,  the  diffusion  assembly  consisted  of  three  distinct  sections:  the  influent 
chamber,  the  soil  column,  and  the  effluent  chamber.  All  components  of  the  apparatus  were 
constructed  from  aluminum  or  stainless  steel  to  minimize  any  interactions  with  VOC  vapors.  The 
length  of  the  soil  column  was  variable,  ranging  from  approximately  0.5  to  24  inches  (1.25  to  61 
cm).  In  the  diffusion  experiments,  the  moisture  content  of  the  aquifer  sand  was  typically  2. 15  to 
2.6  percent  (dry  mass),  and  the  column  was  packed  to  achieve  total  porosities  ranging  from  0.481 
to  0.513. 


Figure  2.  Diffusion  Column  Assembly:  Elevation  View. 


A  cross  section  of  the  diffusion  apparatus  influent  chamber  is  detmled  in  Figure  3.  VOCs  in 
aqueous  solution  were  pumped  at  a  steady  rate  into  the  chamber  at  point  B.  The  drain  assembly 
(H)  provided  both  a  steady  liquid  level  within  the  cell  and  a  “gas  seal”  preventing  any  outflow  of 
vapor.  In  the  gas  headspace  within  the  chamber,  contaminant  vapors  accumulated  as  mass  transfer 
of  VOCs  from  dissolved  phase  to  gas  phase  occurred.  The  VOC  vapor  concentrations  attained 
below  the  aluminum  and  steel  soil  support  screens  (E  and  F)  provided  the  driving  force  for 
molecular  diffusion  through  the  soil.  The  magnitudes  of  the  driving  concentrations,  or  C^,  were 
manipulated  by  altering  either  the  aqueous  solution  concentration  or  by  changing  the  liquid  input 
flow  rate.  Influent  vapor  concentrations  were  measured  during  experiments  by  extracting  0.5  to 
1.0  mL  gas  samples  at  the  sampling  port  (C),  which  was  sealed  with  a  Teflon  backed  silicon 
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septum  and  aluminum  crimp  cap.  The  samples  were  taken  using  Pressure-Lok'^''^  syringes  from 
Precision  Sampling  Corporation.  Samples  were  analyzed  with  a  Hewlett-Packard  5890  Series  II 
gas  chromatograph  that  was  equipped  with  a  stainless  steel  column  (10  ft.  x  1/8”)  packed  with 
20%  SP-2100'r*^  and  0.1%  Carbowax-1500  on  100/200  mesh  Supelcoport™  (Supelco,  INC.; 
Bellafonte,  PA),  a  flame  ionization  detector  (FID)  and  connected  to  a  Hewlett-Packard  3396  Series 
II  integrator.  Influent  vapor  concentrations  ranged  from  approximately  1.5  to  20  percent  of 


Figure  3.  Diffusion  Column  Assembly:  Influent  Chamber,  Elevation  View. 


Figure  4  shows  both  plan  and  elevation  views  of  the  exit  or  effluent  chamber.  This  component  of 
the  diffusion  assembly  provided  a  means  for  maintaining  low  VOC  concentrations,  relative  to 
influent  concentrations;  at  the  downstream  end  of  the  soil  column.  A  thin  stream  of  air  at  100 
percent  relative  humidity  was  directed  through  the  sweep  gas  passage  (B)  at  a  nominal  flow  rate  of 
100  mL/min.  This  "sweep  flow"  carried  contaminants  emerging  from  the  upper  face  of  the  soil 
column  out  of  the  cell  through  the  gas  outlet  at  point  D.  VOC  concentrations  in  the  effluent  sweep 
gas  were  monitored  over  time  via  syringe  extraction  of  0.5  mL  to  1  mL  of  effluent  gas  at  a 
sampling  port  downstream  of  the  chamber.  These  time-varying  effluent  concentration  profiles,  or 
breakthrough  curves,  were  used  to  assess  transient  transport  behavior  in  a  number  of  different 
experiments. 
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Plan  View 


0.10" 
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0.75" 
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\  ^  2.0"(i08cm)  :  . 

;  6.0"  (15.24  cm)  I 

Elevation  View 

Figure  4.  Diffusion  Column  Assembly:  Effluent  Chamber,  Plan  and  Elevation  Views, 


A  schematic  diagram  of  the  entire  diffusion  apparatus  appears  in  Figure  5.  With  the  exception  of 
the  compressed  air  source  tank  (A),  all  components  of  the  experiment  were  maintained  at  25  ± 
0.2  *C  within  a  constant  temperature  chamber. 
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=  Gas  flow  path 
=  Liquid  flow  path 


A.  Compressed  air  and  regulator 

B.  Flow  meter  and  needle  valve 

C.  Hydrocarbon  trap 

D.  Gas  washing  bottles  (water) 

E.  Diffusion  cell  assembly 

F.  Gas  sampling  tee 

G.  Wet-test  flow  meter 

H.  Solvent  source  reservoir 

I.  HPLC  solvent  pump 

J.  Liquid  waste  container 


Figure  5.  Schematic  of  Diffusion  Test  Cell  Assembly  (after  Peterson  et  al.  1988). 


At  steady  state,  the  diffusive  flux  of 
described  by  Fick's  first  law: 


a  compound  through  a  soil  column  of  length  L  (cm)  is 


jC-Co) 

L 


where  J  is  the  vapor  flux  (g/cm2/min),  C  is  the  concentration  of  the  compound  in  the  effluent 
sweep  gas  (g/mL),  Co  is  the  concentration  of  the  compound  in  the  source  vapor  (g/mL),  and  Deff  is 
the  effective  diffusion  coefficient  for  the  VOC  under  the  experimental  conditions  (cm^/min). 


At  steady  state,  the  vapor  flux  through  the  soil  column  may  also  be  determined  via  mass  balance  at 
the  upper  surface  of  the  soil..  The  rate  of  contaminant  mass  flux  through  the  soil  surface,  J,  is 
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equal  to  the  rate  of  mass  exiting  the  diffusion  cell  with  the  sweep  gas  flow; 


where  Q  is  the  sweep  gas  flow  rate  (mL/min),  C  is  the  concentration  of  the  compound  in  the 
effluent  sweep  gas  (g/mL),  and  A  the  cross-sectional  area  of  the  soil  column  (cm2). 


Setting  the  above  two  equations  equal  to  one  another  and  solving  for  Deff  gives: 


QCL 

A(C-Co)’ 


In  analyzing  the  results  of  the  diffusion  experiments,  values  of  source  concentration,  effluent 
sweep  gas  concentration,  and  sweep  gas  flow  rate  data  were  used  to  compute  values  of  Deff.  Data 
points  from  times  near  the  end  of  each  experimental  run  were  selected  for  these  calculations  based 
on  evaluation  of  the  data  and  allowed  adequate  time  for  the  diffusion  process  to  reach  steady  state. 


Advection  experiments 

As  the  primary  means  for  investigating  VOC  mass  transfer,  advective  transport  experiments  were 
conducted  using.  TCE  in  a  bench-scale  column  reactor.  Experimental  details  are  described  by 
Lorden  (1996)  and  are  briefly  summarized  here.  All  experiments  were  conducted  at  25  °C.  A 
stainless  steel  column  (2.36  cm.  I.D.  x  38  cm.  long)  was  filled  with  the  aquifer  sand  and  allowed 
to  equilibrate  to  constant  weight  with  with  a  flowing  gas  stream  at  90%  relative  humidity.  Glass 
saturators  containing  either  water  or  liquid  TCE  were  used  to  prepare  the  gas  stream  to  be  used  in 
experiments.  By  adjusting  the  temperatures  of  the  saturators  and  the  gas  flow  through  each,  it  was 
possible  to  control  the  relative  humidity,  vapor  phase  VOC  concentration,  and  total  gas  flow 
upstream,  of  the  column.  Experiments  were  conducted  at  two  relative  pressures  of  TCE  [10%  P/Pg 
(52.5  mg/L),  and  90%  P/Pg  (472.5  mg/L),  designated  LC  and  HC  for  low  and  high  concentration, 
respectively]  and  two  different  gas  flow  rates  (40  and  80  mL/min,  designated  LF  and  HF  for  low 
and  high  flow,  respectively).  Of  the  two  TCE  relative  vapor  pressures  used,  90%  P/Pg  is 
considered  representative  of  the  “near-field”,  close  to  a  nonaqueous  liquid  phase  VOC  source, 
while  10%  P/Pg  is  more  representative  of  the  “far-field”,  at  contaminated  sites.  The  gas  flow  rates 
were  chosen  so  that  the  resulting  pore  velocity  of  the  gas  would  fall  within  the  range  of  0.001  m/s 
to  0.008  m/s.  This  range  is  thought  to  be  representative  of  the  range  of  velocities  attained  during 
the  soil  vapor  extraction  remediation  process  (Grathwohl  and  Reinhard,  1993).  A  schematic  of  the 
experimental  apparatus  for  the  column  experiments  is  given  in  Figure  6.  In  the  first  section  of  the 
apparatus,  the  gas  flow  rate,  relative  humidity,  and  partial  pressure  of  TCE  were  established.  The 
middle  section  contained  sand-filled  column  through  which  the  gas  flow  traveled.  Provisions  were 
made  to  measure  the  pressure  drop  across  the  column  and  the  relative  humidity  of  the  gas.  Figure 
7  shows: a  detailed  schematic  of  the  column  end  fittings.  The  third  section  of  the  apparatus 
consisted  of  a  gas  chromatograph  with  an  auto  sampling  valve  that  was  used  to  measure  effluent 
TCE  concentrations. 
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Section  1 

1.  Compressed  Air  Cylinder 

2.  Tracer  Gas  Cylinder 

3.  Mass  Flow  Controllers 

4.  Insulated  Water  Baths 

5.  Immersion  Heaters 

6.  SahiratoTS  (TCE) 


=  Allowable  Gas  Flow  Path 
=  Allowable  Water  Flow  Path 


_ I  Section  2 

7.  Saturators  (Water) 

8,  Water  Circulator 

9-  Temp-  Equilibration  Loop 

10.  Column  of  Sand 

11.  Column  Bypass  Line 

12.  Relative  Humidity  Monitor 


*  Section  3 

13.  Gas  Chromatograph 
14-  Sample  Loop  and  Switching 
Valve 

15.  Chromatography  Column 

16.  Flame  Ionization  Detector 

17.  Wet  Test  Meter 


Figure.  6.  Schematic  Diagram  of  Apparatus  for  Study  of  Advective  Transport  of  TCE. 


1-Inch  Cap 


Figure  7.  Schematic  of  Column  Configuration  at  Inlet  and  Outlet. 
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Once  desired  experimental  conditions  had  been  established,  a  breakthrough  profile  (BP)  was 
generated  by  introducing  gas  flow  to  the  column.  After  the  relative  TCE  concentration  (C/Cq)  in 
the  column  effluent  reached  0.99,  the  column  was  allowed  to  equilibrate  with  the  flowing  gas  for 
an  additional  7  days,  during  which  conditions  were  periodically  monitored  to  verify  stability.  The 
column  was  then  isolated  from  gas  flow  for  1  to  2  hours  while  TCE  concentrations  upstream  of  the 
column  were  reduced  to  negligible  levels.  The  gas  flow  was  then  re-directed  through  the  column 
and  a  TCE  elution  profile  (EP)  was  generated.  Effluent  concentrations  were  monitored  until  C/Cq 
was  less  than  or  equal  to  0.001. 

Mathematical  modeling  of  the  data  was  performed  by  Dr.  W.  Chen  using  a  transpon  model  that 
invoked  a  Gamma  distribution  of  sorptive  rate  constants.  The  reader  is  referred  Chen  and  Wagenet 
(1995)  for  a  derivation  and  description  of  the  model  governing  equations.  Methane  breakthrough 
profiles  were  used  to  determine  porosity,  gas  pore  velocity,  and  mechanical  dispersion  for 
subsequent  use  in  mathematical  modeling  of  the  TCE  experiments.  The  TCE  brealrthrough  profile 
at  low  flow  and  low  TCE  relative  pressure  (LC-LF-BP)  was  fit  to  obtain  the  parameters  that 
described  the  Gamma  distribution  of  rate  constants,  a  and  p.  Predictions  were  then  made  of  all 
remaining  TCE  breakthrough  profiles  (BP)  and  elution  profiles  (EP)  at  both  low  and  high  gas  flow 
with  no  adjustable  parameters.  The  TCE  sorptive  distribution  coefficient  used  in  the  model 
computations  was  calculated  directly  from  the  experimentally  measured  TCE  adsorption  isotherm. 

To  accommodate  multiple  sofbates,  the  advective  flow  apparatus  was  modified  as  shown  in  Figure 
8.  The  modifications  permitted  the  presence  of  a  second  vapor  (CF),  increased  tem.perature  (and 
consequently  concentration)  control,  and  computerized  data  acquisition. 

Two  types  of  experiments  were  conducted  and  compared:  single  sorbate  experiments  with  TCE, 
designated  by  “T,”  and  binary  experiments  with  CF  and  TCE,  designated  by  “CT"’.  The 
experimental  designations  below  also  include  a  number  to  indicate  the  relative  pressure  of  the 
sorbates.  Higher  numbers  signify  higher  concentrations  of  TCE  in  the  single  sorbate  experiments 
and  of  CF  in  the  binary  experirhents,  in  which  TCE  concentrations  were  held  constant.  An 
exception  was  experiment  CT5,  which  replicated  experiment  CT4. 

All  concentration  data  was  collected  using  a  gas  chromatograph  (GC)  equipped  with  a  air-actuated 
sampling  valve.  A  computer  program  designed  for  this  research  allowed  full  control  over  sampling 
interval  and  peak  area  integration  parameters  (Weber-Shirk,  1996).  The  experimental  protocol,  in 
general,  was  as  follows.  First,  the  desired  VOC  vapor  concentrations  at  90%  relative  humidity 
were  established  by  varying  both  the  temperatures  and  the  carrier  gas  flow  of  the  gas  washing 
bottles  containing  reservoirs  .of  TCE,  CF,  and  water,  respectively.  The  column  of  aquifer  sand 
was  then  exposed  to  the  vapors  for  two  hours.  VOC  breakthrough  typically  occurred  within  the 
first  10  minutes.  At  the  end  of  the  two  hour  equilibration  the  column  was  isolated  and  the 
plumbing  purged  to  remove  any  residual  VOCs.  Clean,  humidified  air  was  then  used  to  elute  the 
column  in  three  stages,  separated  by  intervals  18  and  6  hours,  respectively.  In  each  stage  the 
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effluent  peak  areas  were  monitored  until  the  gas  concentrations  approached  the  GC  peak  area 
detection  limit.  The  interval  between  elutions  allowed  VOC  concentrations  in  the  column  to 
rebound  by  slow  desorption  from  the  stationary  phase.  Tests  indicated  that  this  procedure 
permitted  the  recovery  of  more  than  99%  of  the  VOC  mass  sorbed  onto  the  column. 


1.  Source  Ai  r  or  Tracer  Gas  Cylinder 

2.  Dessicator 

3.  Mass  Flow  Controller  Transducer 

4.  Mass  Flow  Controller 

5.  Constant  Temperature 

6.  VOC  Gas  Washing 

7.  WaterGasWashiI:^g 


8.  Cor\st ant  Temperature  Circulator 

9.  Soil  Column 

10.  Manometer 

11.  Flowmeter  with  Restrictor 

12.  Gas  Chromatograph  (GC) 

13.  Computer 

14.  PurgingCylinder  (N,  Gas) 


Figure  8.  Schematic  of  the  Modified  Advective  Flow  Apparatus. 
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Results  and  Discussion: 

VOC  Sorption 

At  low  relative  pressures,  TCE  adsorption  was  described  by  a  linear  isotherm  (see  figure  9). 


Figure  9.  Linear  Isotherm  for  Sorption  of  TCE  Vapor  on  Moist  (M.C.  ==  A%)  Aquifer 
Sand. 


The  soil-vapor  distribution  coefficient,  ,  was  determined  from  the  following  equation  derived 
from  mass  balances  by  Ong  and  Lion  (1991); 


where  C  denotes  concentration,  V  denotes  volume,  and  the  subscripts  B  and  S  denote  control 
(blank)  and  sample  vials,  respectively.  M  is  the  mass  of  a  moist  soil  sample  minus  the  mass  of 
associated  water. 


The  value  of  K^'  is  given  by  the  slope  of  a  plot  of  VgC^/V^Cs  versus  M/l^  •  es  shown  in  Figure  9 
and  equaled  0.60  ±  0. 1 1  ml/g  at  the  95%  confidence  interval. 

Ong  and  Lion  (1991)  developed  the  following  predictive  equation  for  K^'  which  is  applicable 
when  the  average  number  of  bound  molecular  layers  of  water  exceeds  five,  the  soi;  water  behaves 
ideally,  and  when  VOC  condensation  and  sorption  at  the  gas-water  interface  are  negligible: 

a:;  =  aTj  /  /:„  +  (m.c.,  %)/( i  oo  ) 
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where  JTj,  is  the  Henry’s  constant  for  the  compound,  is  the  saturated  distribution  coefficient 
for  sorption  from  the  aqueous  to  the  solid  phase  and  M.C.  is  the  sorbent’s  moisture  content. 

Substituting  =  0.20  ml/g  (Lorden,  1996),  M.C.  =  4.0%,  and  K.^  =  0.397  (Garbarini  and 
Lion,  1985)  yields  -  0.60  ml/g.  At  this  level  of  significant  figures,  the  predicted  value  is  the 
same  as  the  experimentally  determined  value.  The  excellent  agreement  between  experiment  and 
theory  suggests  that  the  value  of  detennined  in  this  research  is  accurate,  and  that  model  of  Ong 
and  Lion  (1991)  can  be  used  to  estimate  for  a  given  moisture  content. 

Non-linear  adsorption  isotherms  were  developed  for  both  water  and  TCE  (Figures  10  and  11) 
using  the  gravimetric  apparatus  described  by  Ong  and  Lion  (1991).  The  adsor|:)f.ion  of  water  on 
oven-dried  sand  yielded  an  isotherm  with  characteristics  fitting  the  Type  I  BET  isotherm  described 
by  Gregg  and  Sing  (1982).  . 
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Figure  10.  Water  Sorption  on  Oven-Dried  Aquifer  Material. 
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The  TCE  isotherm  on  the  moist  sand  (Figure  11)  has  characteristics  comparable  to  those  of  the 
Type  V  isotherm  as  defined  by  Gregg  and  Sing  (1982).  The  aquifer  sand  sorbent  had  additional 
capacity  for  TCE  even  at  a  relative  humidity  of  90%,  suggesting  that  condensation  of  TCE  could 
occur  even  in  nearly  saturated  soils.  The  TCE  adsorption  isotherm  also  indicate?  that  condensed 
vapor  can  constitute  a  high  fraction  of  the  TCE  sorbed  at  high  relative  pressures  of  TCE. 
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Figure  11.  TCE  Sorption  on  Aquifer  Material  at  90%  Relative  Humidity. 

Competitive  sorption  between  VOCs  at  high  relative  vapor  pressures  was  also  observed.  For 
example,  this  was  evident  in  the  case  of  PCE  in  the  presence  of  CF.  As  illustrated  in  Figure  12  the 
PCE  isotherm  was  non-linear  and  PCE  sorption  was  suppressed  in  the  presence  of  CF. 


Figure  12.  Sorption  Isotherms  for  PCE  on  Moist  Aqufer  Material  in  the  Presence  and 
Absence  of  CF. 
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Advective  Transport 

The  breakthrough  curves  shown  in  Figure  13  demonstrate  the  reproducibility  of  replicate  column 
experiments. 
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Figure  13.  Replicate  Breakthrough  Profiles  for  TCE  P/Pq  =  0.90,  Q  =  40.0  mL/min 


Previous  research  has  suggested  that  a  Gamma  distribution  of  desorption  rates  can  serve  as  a 
robust  tool  for  characterization  of  contaminant  mass  transfer  (Connaughtcn  et  gL,  1994;  Chen, 
1994;  Culver,  et  al.  1997).  Gamma  model  parameters  of  a  =  1.85  and  [3  ==  2.22  x  10^  minutes 
were  obtained  from  a  fit  of  the  LC-LF-BP  column  experiment.  The  value  of  the  distribution 
coefficient  used  for  the  fit  was  11.8  ml/g,  and  was  calculated  directly  from  the  TCE  adsorption 
isotherm  at  a  relative  pressure  of  10%  P/Pq  (see  Figure  1 1).  The  fitted  values  obtrined  for  a  and  P 
suggest  that  slow  sorption  sites  predominated  in  the  advective  transport  expenrnents.  The  mean 
first  order  rate  transfer  coefficient,  k  ,  was  equal  to  8.33  x  10’5  min'b  A  prcb.ability  density 
function  generated  from  the  fitted  values  of  a  and  p  is  shown  in  Figure  14.  Note  that  a  peak  in  the 
distribution  occurs  at  k  4.0  x  10-5  min'l,  and  that  the  mean  k  value  lies  sligh  tly  to  the  right  of  this 
value. 
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Figure  14.  Experimentally  Determined  Gamma  Distribution  of  First  Order 
Applicable  to  Sorption  and  Desorption  of  TCE  over  a  Wide  Ra 
Pressures  and  at  Two  Gas  Flow  Rates. 


Figure  15  illustrates  the  the  predictions  of  the  Gamma  site  model  for  the  breal'tthro 
of  the  LC-HF  experiment.  Figure  16  provides  a  comparison  for  the  prediction; 
model  to  the  elution  profile  (EP)  from  the  LC-HF  experiment.  The  agreement  bet 
and  the  observations  is  excellent.  Model  predictions  were  also  obtained  with  the  s 
and  P  for  all  other  experimental  conditions  including  those  at  a  high  gas  concentrat 
flow  velocity  and  agreed  well  with  the  data  [see  Lorden  (1996)  for  additions’  cor 
model  and  experimental  data].  Additional  calculations  (in  progress)  suggest  tlrat  th 
predictions  are  significantly  superior  to  those  of  two  other  commonly  invoked  par 
site/two  region  model  and  a  pore  diffusion  model. 
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Figure  15.  Comparison  of  Experiirieiital  Data  with  the  Predictions  of  the  Gsmrr  a  Mode!  for  the 
Breakthrough  Profile  of  the  Low  Concentration,  High  Fiovv  Coh-inm  Experiment. 
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Figure  16. 


Comparison  of  Experimental  Data  with  the  Predictions  of  the  Gamr!  a  Model  for  the 
Elution  Profile  of  the  Low  Concentration,  High  Flow  Column  E  xperiment. 
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A  listing  of  the  experimental  conditions  from  column  experiments  in  which  the  behavior  of 
mixtures  was  considered  is  given  in  in  Table  1. 


Table  1:  Experimental  Conditions  for  Comparison  of  TCE  and  TCE  with  GF  Elution  Profiles 


Experiment 

Relative  Vapor  Pressure 

Mass  So- 

-bed 

(  mg! kg  zoei 

sand) 

TCE 

CF 

TCE 

CF 

T1 

0.10 

0.0574 

T2 

0.10 

- 

0.0616 

- 

T3 

0.10 

- 

0.0557 

- 

T4 

0.23 

- 

0.1037 

- 

T5 

0.33 

- 

0.1815 

- 

CTl 

0.10 

0.03 

0.0666 

0.C413 

CT2 

0.10 

0.06 

0.0518 

0.C693 

CT3 

0.10 

0.10 

0.0527 

0.1299 

CT4 

0.10 

0.21 

0.0545 

0.2635 

CT5 

0.10 

0.21 

0.0531 

0.2500 

CT6 

0.10 

0.37 

0.0510 

0.4  5S0 

Based  on  the  observations  from  the  column  studies  listed  in  Table  1,  pseudo  isotherms  for  TCE 
and  CF  may  be  developed  as  shown  in  Figures  17  and  18,  respectively.  These  isotherms  display 
the  amount  of  vapor  uptake  on  the  column  for  a  given  influent  VOC  concent:  ation  and  are  subject 
to  possible  mass  transfer  limitations  that  may  have  restricted  vapor  contac  with  the  stationary 
phase  within  the  2  hour  equilibration  period  that  was  used. 


The  sorption  data  for  both  compounds  in  the  range  of  concentrations  examined  were  best  fit  to 
linear  isotherms.  However,  the  data  may  also  be  fit  to  BET  type-IIT  isotherms,  winch  are  linear  at 
lower  vapor  concentrations,  with  virtually  identical  R2  values.  This  obsers  ation  is  consistent  with 
the  results  obtained  by  other  researchers  for  VOC  sorption  under  conditions  of  ad  'ective  flow  and 
90%  relative  humidity  (Pennell  et  al,  1992;  Thibaud  et  al,  1993). 

The  TCE  vapor  concentration  was  held  constant  at  a  relative  vapor  pres.sure  of  0. 10  at  25  ^C 

iii  experimental  runs  T1-T3  and  all  the  mixture  experiments  (CT  runs  1-6:.  Tf  e  mean  of  TCE 
sorbed  in  the  single  sorbate  experiments  was  0.0582  mg  TCE/kg  et  sand,  compared  to  0.0526 
mg/kg  in  the  binary  runs  (see  Table  2).  [Note,  the  mass  of  TCE  sorbed  in  experiment  CTl  is 
considered  ah  anomaly  due  to  a  temperature-induced  problem  with  the  TCE  vape  r  Eow.]  A  non- 
parametric  statistical  analysis,  the  Wilcoxon  rank  sum  test,  was  applied  to  the  TCE  sorption  data. 
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.  The  results  indicate  that  the  average  mass  of  TCE  sorbed  is  signifiean:Iy  d  ffercnt  in  the  binary 
experiments  to  the  95%  confidence  level. 


Partial  Pressure  (P/Pp) 


Figure  17.  TCE  Isotherms  Based  on  Column  Experiments. 


Figure  18.  CF  Isotherms  Based  on  Column  Experiments, 
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The  results  indicate  that  the  average  mass  of  TCE  sorbed  is  significantly  different  in  the  binary 
experiments  to  the  95%  confidence  level. 


Figure  17.  TCE  Isotherms  Based  on  Column  Experiments. 


Figure  18.  CF  Isotherms  Based  on  Column  Experiments. 
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The  effect  of  C3^  on  TCE  sorption  was  not  dependent  on  the  CF  vapor  concentration,  but  rather  a 
presence/absence  relationship.  CF  concentrations  spanning  an  order  of  magnitude  reduced  TCE 
sorption  by  approximately  the  same  amount,  an  average  of  0.0056  mg/kg.  Th is  finite,  quantifiable 
effect  suggests  that  CF  completely  may  effectively  compete  with  TCE  fo.  a  particular  class  of 
sorption  site(s)/ 

There  arc  five  types  of  sorption  “sites”  (or  uptake  mechanisms)  that  are  gene-ally  considered 
applicable  to  vapor-phase  sorpdcm  in  the  vadose  zone:  1)  VOC  dissolutio.  in  sci!  moisture;  2) 
VOC  partitioning  from  soil  moisture  into  soil  organic  matter;  3)  vOC  c  ’sorption  at  mineral 
surfaces;  4)  VOC  diffusion  into  the  pores  of  soil  grains;  and  5)  VOC  son'. ion  at  the  gas-water 
interface.  Both  VOC  dissolution  and  partitioning  arc  typically  linear  prccessc  s  an  i  are  not  subject 
to  con4)etitive  interactions  between  VOCs.  Since  water  effectively  out-comperes  not  ionic  organics 
such  as  TCE  and  CF  for  polar  mineral  surface  sites,  any  competitive  VOC  sorption  at  the  mineral 
surface  would  be  expected  to  be  of  a  negligible  magnitude.  Migration  of  VOCs  into  intraparticle 
pores  is  governed  by  aqueous  diffusion,  a  slow  transport  mechanism.  Given  the  experimental 
equilibration  time  of  two  hours,  coupled  with  tihe  fact  that  over  99%  mass  recover/  of  sorbed  YCXi: 
was  accomplished  in  three  elution  phases  spanning  24  hours,  it  seems  extremely  unlikely  that  a 
significant  annount  of  VCX^  mass  was  located  in  intraparticle  pores.  Therefo  e,  it  s  deemed  likely 
that  the  observed  competition  between  CF  rad  TCE  must  have  occurred  at  the  gas-  vater  interface. 

ScHne  researchers  have  det^mined  that  sorption  at  the  gas-waier  interface  may  be  signiHcant.  and 
perhaps  as  much  as  50%  of  the  total  under  certain  conditions  (Pennell  et  oi,  IS  92).  Karger,  et 
ai  (1971a)  found  vaptM’-phase,  interfacial  adsorption  to  be  important  under  advective  conditions. 
In  another  report,  Karger  rad  colleagues  (1971b)  suggested  that  it  was  possible  for  compounds 
^rbed  at  the  gas-water  interface  to  fonn  hydrogen  bonds,  and  that  interfacia  adscrp  .ion  increased 
with  increasing  polarizability  of  the  adsorbing  species.  Valsaraj  (1994)  a.  vance d  this  concept 
further  by  demonstrating  that  adsorption  at  the  interface  was  thermodyr  a  >  icaliy  favorable  for 
hydn^obic  compounds. 


Given  the  similar  experimental  conditions  in  this  research  to  those  of  the  other  iir.  es 
above  and  the  hydrophobicity  of  TCE  rad  CF,  it  seems  likely  that  significa  sor  •*; 
water  interface  occurred.  Corley  (1996)  found  that  the  fraction  of  the  ior<  1  ;ri.ass  sc 
attributable  to  adsorption  at  the  gas-water  interface  declined  with  increasing  vat:  c 
the  ratio  of  interfacial  area  to  volume  (of  water)  decreases  at  higher  moi  ture  cc 
dissolution  into  soil  moisture  becomes  increasingly  important.  Hoff  sr  al.  ' ;  993' d  ■ 
the  fractional  mass  adsorbed  at  the  gas-water  interface  declined  with  increa ?i  g  ac  cil 
This  is  not  surprising  given  that  sorption  into  intraparticle  pores  is  diffu  ion-lici 
becomes  increasingly  inqKjrtant  at  longer  equilibration  times.  While  neither  c  f  the  e  : 
definitively  prove  that  tire  gas-water  interface  has  a  relatively  small,  finite  n  n  ,ber  f 
they  do  suggest  that  this  is  likely.  If  this  is  true,  then  competiti  on  het v.  -en  C  F 
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interfacial  sorption  sites  seems  probable.  Since  CF  has  both  greater  polarizability  and  has  twice  the 
hydrogen  bond'  electron  donor  basicity  as  TCE,  it  should  be  an  effective  competitor  for  interfacial 
,  sorption  sites,  as  the  data  from  this  research  suggests. 

In  order  to  examine  competitive  effects  on  the  relative  removal  rates  of  the  VOCs,  the  mass 
removed  from  the  soil  column  over  time  during  the  first  elation  phases  of  each  experiment  was 
normalized  by  the  total  mass  removed  over  all  three  elutions.  This  allowed  for  column  experiments 
with  different  total  masses  of  VOCs  sorbed  to  be  compared  on  a  the  basis  of  percent  mass 
remaining  in  the  column.  The  transformed  data  was  plotted  against  the  number  of  pore  volumes  of 
clean  air  passed  through  the  column.  The  removal  of  TCE  in  experiments  at  which  TCE  was  held 
at  a  constant  vapor  concentration  of  P/Pq  =  0.10  is  shown  in  Figure  19. 


Figure  19.  TCE  elution  in  experiments  with  TCE  P/Pq  =  0.10. 

Figure  19  suggests  that  the  sorbed  TCE  sorbed  was  removed  faster  in  the  presence  of  CF,  and  that 
increasing  the  amount  of  CF  further  increased  the  removal  rau'.  Since  the  number  of  moles  of  TCE 
sorbed  was  roughly  constant  ih  each  case,  there  also  appears  to  be  a  tread  that  indicates  that  the 
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removal  rate: of  TCE  was  dependent  on  the  total  number  of  moles  ofboth  VOCs  sorbed  (see  Table 
2).  This  trend  is  also  seen  in  die  data  for  the  removal  of  CF  in  the  bmar>'  experimenLs.  As  shown 
in  Figure  20,  CF  removal  increased,  on  a  percentage  basis,  with  an  increasing  number  of  total 
moles  of  VOC  sorbed.  This  same  trend  is  seen  in  the  single  sorbate  experiments  involving  only 
TCE  (Figure  21).  Thus  the  effect  of  the  presence  of  CF  on  the  desorption  removal  rate  of  TCE 
can  not  be  discerned  from  the  apparent  effect  of  the  increase  m  the  total  number  of  moles  of  VOCs 

sorbed  to  the  sand. 


Figure  20,  CF  Removal  in  Column  Experiments. 
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Table  2>  Total  Number  of  Moles  Sorbed  in  Each  Experiment. 


Experiment 

Partial  Pressure 

TCE  CF 

Moles  Sorbed 
(mmol) 

TCE  CF 

Total  Moles  Sorbed 
(mmol) 

T1 

■i 

. 

2.190 

2.190 

T2 

- 

2.350 

- 

2.350 

T3 

BB 

- 

2.124 

- 

2.124 

T4 

- 

3.954 

- 

3.954 

T5 

■H 

- 

6.922 

- 

6.922 

CTl 

0.10 

0.03 

2.538 

1.431 

3.969 

CT2 

0.10 

0.06 

1.976 

2.402 

4.378 

CT3 

0.10 

0.10 

2.009 

4.501 

6.509 

CT4 

0.10 

0.21 

2.084 

9.131 

11.215 

CT5 

0.10 

0.21 

2.025 

8.663 

10.688 

CT6 

0.10 

i  0.37 

1.944 

15.872 

17.816 

JUL  lite;  'y/  wiiJibrn 


r.c: 


Diffusive  Transport  of  VOCs 

Diffiision  experiments  involving  a  number  of  different  VOC  concentrations  and  combinations  were 
conducted.  For  convenience,  each  experiment  is  described  with  reference  to  maximum  influent 
concentratidns  measured  for  each  compound.  A  brief  summary  of  concentration  conditions  is 
shown  in  Table  3.  In  the  table,  and  subsequently  in  this  report,  experiments  will  be  referred  to  by 
listing  the  first  letter  of  each  compound  involved,  in  alphabetical  order,  followed  by  the  run 
number.  For  example,  the  sixth  experiment  involving  TCE  and  CF  is  designated  CT-6. 

An  example  of  raw  data  collected  during  experiment  CT-7  is  displayed  in  Figures  22  and  23.  The 
CF  and  TCE  influent  profiles  shown  in  Figure  22  display  a  gradual  approach  to  a  maximum 
concentration  which  was  a  function  of  the  influent  liquid  flow  rate  and  the  mass  transfer  of  VCKls 
from  the  liquid  phase  to  the  gas  phase  in  the  influent  chamber.  Corresponding  effluent  profiles  in 
Figure  23  show  a  1  to  2  hour  lag  between  the  initiation  of  the  experiment  and  the  appearance  of 
VOCs  in  the  effluent  sweep  gas.  This  delay  is  consistent  with  retardation  of  the  diffusive  vapor 
flux  by  the  aquifer  sand  medium. 


Table  3.  Diffusion  Experiment  Concentration  Summary 


Experiment 

CF  Influent 
Concentration 
:(%  Of  Psat) 

PCE  Influent 
Concentration 

(%  of  Psat) 

TCE  Influent 
Concentration 
(%  of  Psat) 

C-3 

2.37 

-- 

— 

C-4 

21.05 

— 

— 

CT-5 

2.20 

— 

1.67 

CT-6 

19.93 

- 

1.44 

CT-7 

18.64 

— 

3.28 

CT-8 

16.04 

— 

13.56 

PT-5 

— 

2.78  ‘ 

1.36 

T-15 

— 

3.78 

T-16 

— 

1.54 

T-18 

— 

— 

1.46 

T-19 

— 

— 

13.00 
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Figure  22.  Experiment  CT-T  Influent  Concentration  Profiles. 
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,A  number  of  studies  have  been  conducted  to  determine  the  relationship  between  physical 
cha,racteristick  of  porous  media  and  vapor  diffusion  rates.  One  widely  accepted  model  which 
accounts  for  physical  impedance  to  diffusion  was  proposed  by  Millington  and  Quirk  (1961).  The 
model,  based  on  the  concept  of  porosity  and  the  probability  of  connected  pores,  provides  a 
relationship  between  vapor  diffusion  coefficients  in  air  and  in  a  porous  media.  The  form  of  the 
Millington-Quirk  model  is: 


Were,  Dq  is  the  free  air  diffusion  coefficient  (cm^/min),  a  is  air-filled  porosity,  <{>  is  total  porosity, 
and  Dgff  is  the  "effective",  dr  physically  hindered,  diffusion  coefficient.  Values  of  D^ff  are 
independent  of  sorption  or  any  other  reaction. 

..  As  described  above  under  meftods,  the  effective  diffusion  coefficient  for  a  vapor  can  be  calculated 
from  steady  state  diffusion  data  using  Pick’s  first  law.  Concentration  gradients  observed  during 
the  late  stages  of  several  diffusion  experiments  (after  influent  and  effluent  concentrations  stabilized) 
were  used  to  compute  Dgff  values.  These  values  were  compared  to  predictions  made  using  the 
Millington-Quirk  model.  The  experimentally  derived  and  empirically  predicted  Deff  values  showed 
very  close  agreement. 

Since  the  Millington-Quirk  model  was  shown  to  serve  as  an  accurate  means  for  describing  the 
physical  hindrances  to  VOC  diffusion,  it  was  possible  to  look  more  closely  at  sorption  phenomena. 
This  was  accomplished  using  a  numerical  simulation  to  predict  the  transient  behavior  or  vapor 
breakthrough  in  the  diffusion  experiments.  For  the  purpose  of  this  report,  a  very  brief  description 
of  the  model  is  presented.  Details  pertaining  to  the  numerical  model,  along  with  a  complete 
analysis  of  model  output,  may  be  found  in  Gilbert  (1997,  M.S.  thesis,  in  preparation). 


In  general,  the  transient  diffusive  transport  of  a  vapor  through  porous  media  is  described  by  Pick’ 
second  law: 

dt  J 

where  Daon  is  the  “apparent”  diffusion  coefficient  and  t  is  time. 

opp 


In  contrast  to  the  effective  diffusion  coefficient,  the  apparent  diffusion  coefficient  includes  sorption 
related  retardation  effects.  Pick’s  second  law  is  often  expressed  as. 


dt  R  ) 

where  R  is  a  dimensionless  retardation  factor.  The  value  of  R  is  always  greater  than  or  equal  to  1 , 
with  larger  values  corresponding  to  more  significant  sorption  and  correspondingly  slower 
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diffusion.  ^  investigation  of  VOC  sorptive  behavior  in  this  study  required  an  analysis  of  R. 
Consequently,  a  solution  to  the  above  equation  was  necessary. 

Many  analytical  solutions  to  Picks  second  law  have  been  documented  for  particular  boundary 
I  conditions  (Carlslaw  and  Jaeger,  1959;  Crank,  1975).  The  combination  of  boundary  conditions 
encountered  during  this  study,  however,  made  it  necessary  to  utilize  a  numerical  solution.  This 
was  accomplished  by  developing  a  finite  difference  model  that  simulated  the  diffusion  experiments. 
I  Input  parameters  for  the  model  included  D^ff  (determined  using  the  Millington-Quirk  expression), 
a,  (|>,  Q,  the  length  of  the  soil  column,  and  an  expression  for  the  influent  VOC  vapor 
concentrations  as  a  function  of  time.  Given  these  known  parameters  and  an  estimated  value  for  R, 
the  model  provided  a  prediction,  using  no  adjustable  parameters,  of  the  expected  effluent  VOC 
.  concentration  profile. 

Figure  24  shows  the  measured  influent  TCE  data  from  experiment  T-18,  and  a  cur\'e  used  to 
represent  to  the  data.  A  curve  fitting  application  was  used  to  identify  an  equation  which  best 
represented  experimental  measurements.  In  all  cases,  an  expression  of  the  following  form  was 
used: 

C,  =A[l~exp(-Bf)] 

where  A  and  B  are  adjustable  fitting  parameters.  In  all  cases,  this  curve  fitting  process  yielded  R2 
values  of  0.9^8  or  greater.  Influent  expressions  determined  in  this  way  were  used  as  input  to  the 
diffusion  coihputer  model. 
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The  output  from  the  numerical  simulation  of  the  experiment  T-18  TCE  diffusion  is  shown  along 
with' the  measured  effluent  concentration  profile  in  Figure  25.  The  value  of  the  retardation  factor, 
R,  which  appears  in  Figure  24  is  based  upon  the  assumption  of  a  linear  sorption  isotherm  for  TCE, 


0  200  400  600  800  1000  1200  1400  1600 

Elapsed  Time  (min) 

Figure  25.  Experiment  T-18  effluent  TCE;  Data  and  Model  Prediction. 


Referring  to  Table  3,  it  may  be  seen  that  the  maximum  influent  concentration  for  TCE  during 
experiment  T-18  was  relatively  low,  1,46%  of  Pgaf  Sorption  isotherms  displaying  linearity  al  low 
sorbate  vapor  concentrations  have  been  observed  by  several  researchers  (Chiou  and  Shoup,  1985; 
Ong  et  aL,  1991)  With  the  assumption  of  a  linear  sorption  isotherm,  an  expression  for  the 
retardation  factor  may  be  written  as: 


R  =  l  + 


pK 

a 


where  p  is  the  sorbent  bulk  density  (g/cm^),  a  is  the  air-filled  porosity  and  K^j’  (mL/g)  is  the 
overall  partition  coefficient  describing  equilibrium  between  the  vapor  phase  and  the  sorbed  phase. 
As  noted  above,  a  method  for  predicting  Kj’  from  easily  measured  parameters  has  been 
documented  in  the  literature  (Ong  and  Lion,  1991)  and  was  confirmed  by  this  research. 

■  Kh  IOOKh 

where  is  the  liquid  phase  or  saturated  partition  coefficient,  Kh  is  the  dimensionless  Henry’s 
constant,  and  0  is  the  soil  moisture  content  expressed  as  a  percentage  of  dry  soil  mass.  The  value 
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of  Kd  is  both  compound  and  sorbent  specific,  and  was  measured  for  TGE  and  the  sorbent  used  in 
this  study  (Lx)rden,  1996),  Henry’s  constant  is  readily  available  for  a  large  number  of  compounds. 
Sofbent' moisture  content  was  directly  measured  in  this  research.  Knowledge  of  all  variables  on 
:  the  right-hand  side  of  the  aboye  equation  allowed  for  computation  of  K^’.  This  Kj’  value  was  then 
.  used  to  calculate  an  expected  retardation  factor  for  each  experiment. 

The  retardation  coefficient  used  for  the  model  prediction  in  Figure  26  was  calculated  as  described 
above.  Consequently,  the  numerical  modefhad  no  fitting  parameters,  and  never-the-less  displayed 
excellent  agreement  with  measured  effluent  data.  This  suggests  that  the  numerical  model  employed 
is  a  valid  predictive  tool  for  TCE  at  low  influent  concentrations.  This  procedure  was  employed  to 
model  TCE  data  from  all  experiments.  It  was  observed  that  use  of  independently  determined  R 
values  yielded  excellent  model  predictions  of  breakthrough  behavior  for  TCE  at  low 
concentrations,  but  gave  less  accurate  estimates  of  effluent  profiles  at  high  TCE  concentrations  or 
for  experiments  involving  more  than  one  VOC. 

.  The  diffusion  apparatus  designed  for  this  research  permitted  introduction  of  VOC  mixtures  to  the 
influent  chamber;  and  therefore  examination  of  the  impact  of  potential  competition  among  VOC 
vapors  on  transport  behavior,.  Experiments  were  conducted  using  binary  vapor  mixtures  of  CF 
and  TCE.  Experimental  proceiiures  allowed  for  manipulation  of  the  relative  concentrations  of  each 
compound  in  the  mixture.  By  comparing  observed  VOC  transport  behavior  from  binary 
experiments  with  profiles  from  independent  TCE  or  CF  experiments,  assessments  of 
presehce/absence  or  concentration  dependent  impacts  of  a  co-sorbate  were  possible. 

Plots  such  as  the  one  displayed  in  Figure  26  allow  for  qualitative  comparisons  of  VOC  behavior 
under  varying  mixture  or  concentration  circumstances.  Generally  speaking,  experiments  were 
conducted  using  different  maximum  influent,  and  consequently  effluent,  vapor  concentrations. 
Since  it  was  desirable  to  compare  the  time  dependent  effluent  profiles  from  different  experiments,  it 
was  necessary  to  normalize  the  effluent  breakthrough  graphs  for  each  experiment.  This  was 
accomplished  by  dividing  each  measured  effluent  data  point  by  its  respective  maximum  observed 
influent  concentration.  The  result,  as  shown  in  Figure  26,  was  a  vertical  axis  which  had  been 
“normali2ed’\  allowing  for  direct  visual  comparison  of  breakthrough  trends  from  multiple 
experiments. 

Figure  26  shows  normalized  effluent  CF  profiles  from  experiments  with  varying  amounts  of  TCE 
in  the  influent  vapor.  As  the  graph  legend  indicates,  experiment  CT-8  had  an  influent  TCE 
concentration  approximately  four  times  greater  than  that  in  the  influent  of  the  other  experiments 
shown  in  the  figure.  The  departure  of  the  CT-8  data  from  the  other  curves  suggests  more  rapid 
breakthrough  of  CF  (and  hence  reduced  CF  retardation).  It  is  possible  to  explain  this  breakthrough 
behavior  in  terms  of  competition  between  CF  and  TCE  in  sorption  reactions.  If  CF  sorption  is 
inhibited  by  TCE,  less  retardation  of  CF  would  be  expected,  leading  to  more  rapid  transport  and 
breakthrough  of  CF.  The  breakthrough  curves  in  Figure  25  are  also  consistent  with  what  may  be  a 
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, concentration  dependent  effect,  with  higher  levels  of  TCB  having  a  more  pronounced  impact  on  CF 
behavior. 


Figure  26.  CF  Breakthrough  with  Varying  TCE  Concentration 


Several  relationships  of  the  type  depicted  in  Figure  26  were  observed  upon  normalization  of  the 
experimental  data.  However,  since  the  normalization  procedure  did  not  account  for  small 
differences  in  soil  conditions  (e.g.  porosity,  moisture  content)  or  time  dependent  influent  profiles, 
these  visual  observations  were  not  deemed  to  be  conclusive  proof  of  competitive  interactions 
between  VOCs.  A  more  detailed  examination  of  the  experimental  results  is  still  underway,  the 
results  of  which  will  be  included  in  the  thesis  associated  with  this  ongoing  research  (Gilbert,  1997, 
M;S.  thesis,  in  preparation). 


Conclusions 

The  research  described  in  this  report  is  continuing  under  the  auspices  of  a  DOD  augmentation  grant 
“Augmentation  for  Basic  Research  on  Vapor  Transport  in  Soil”  (AFOSR  project  no:  F49620-94-1- 
0316).  As  a  result,  the  full  analysis  of  the  data  and  its  implications  has  not  been  completed. 
Nevertheless,  several  conclusions  can  be  drawn  at  this  stage  as  follow. 

At  low  relative  vapor  pressure,  the  sorption  of  TCE  on  a  moist  sandy.aquifer  material  displayed  a 
linear  isotherm  that  could  be  predicted  using  the  model  of  Ong  and  Lion  (1991).  At  higher  relative 
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pressures,  the  TCE  sorptioii  isotherm  was  non-linear,  in  accordance  with  BET  theory  and 
indicating  vapor  uptake  occurs  at  the  gas-water  interface, 

•  The  Gamma  model,  which  employs  a  statistical  distribution  of  mass  transfer  rates,  has  now  been 
successfully  employed  to  describe  the  advectivc  transport  of  VOCs  in  unsaturated  porous  media. 
Model  parameters  obtained  from  a  single  column  experiment  and  from  equilibrium  sorption 
measurements  were  capable  of  predicting,  with  no  adjustable  parameters,  the  breakthrough  and 
elution  of  TCE  vapor  over  a  wide  rage  of  relative  pressures  and  gas  flow  velocities. 

Column  experiments  comparing  TCE  vapor  sorption  individually  and  as  a  mixture  with  CF  show 
inild  but  significant  suppression  of  TCE  sorption  in  the  presence  of  CF.  Both  TCE  and  CF  elution 
rates  appeared  to  be  enhanced  when  greater  amounts  of  vapor  were  sorbed.  This  result  is 
consistent  with  the  expectation  that  increases  in  vapor  uptake  may  occur  at  the  gas  water  interface 
and  that  the  release  rate  of  this  material  is  faster  that  that  of  dissolved  VOC  or  VOC  that  is  sorbed  to 
soil  organic  matter.  However,  since  vapor  sorption  was  generally  less  than  that  at  equilibrium 
based  on  gas-water  and  water-sand  distribution  coefficients,  the  occurrence  of  VOC  uptake  at  the 
gas-water  interface  can  not  be  definitively  established  in  the  column  experiments.  To  the  extent 
that  VOC  condensation  occurs  in  the  field,  an  expected  consequence  would  be  that  gas  phase  soil 
clean  up  technologies,  such  as  soil  vapor  extraction,  would  to  be  most  efficient  in  removing  at  high 
levels  of  VOC  contamination,  but  the  mass  transfer  rate  for  removal  would  decline  after  removal  of 
VOC  at  the  gas  water  interface; 

The  Millington-Quirk  formulation  for  calculation  of  the  effect  of  soil  tortuousity  on  the  steady  state 
diffusion  coefficient  of  vapors  has  been  shown  to  give  accurate  predictions.  At  low  relative 
pressures,  the  transient  diffusion  coefficients  of  vapors  show  the  expected  retardation  attributable 
to  linear  sorption  processes.  However,  vapor  diffusion  in  mixtures  appears  to  be  enhanced  in 
some  cases,  as  does  vapor  diffusion  at  high  relative;  pressures.  Both  of  these  observations  suggest 
nonrlinear  isotherms  for  vapor  uptake  are  operative  and  that  vapor  competition  and/or  condensation 
at  the  gas  water  interface  can  act  to  enhance  diffusive  transport.  A  consequence  of  this  observation 
is  that  diffusive  spread  of  vapor  mixtures  in  the  near  field  close  to  source  materials  should  happen 
somewhat  more  rapidly  than  would  be  expected  based  on  models  that  employ  linear  \'OC  phase 
distribution  processes  would  predict. 
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